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Abstract 

In this paper, we consider the problem of multi-relay selection for multi-stream cooperative MIMO systems with M relay 
nodes. Traditionally, relay selection approaches are primarily focused on selecting one relay node to improve the transmission 
reliability given a single-antenna destination node. As such, in the cooperative phase whereby both the source and the selected 
relay nodes transmit to the destination node, it is only feasible to exploit cooperative spatial diversity (for example by means 
of distributed space time coding). For wireless systems with a multi-antenna destination node, in the cooperative phase it is 
possible to opportunistically transmit multiple data streams to the destination node by utilizing multiple relay nodes. Therefore, 
' we propose a low overhead multi-relay selection protocol to support multi-stream cooperative communications. In addition, we 

I derive the asymptotic performance results at high SNR for the proposed scheme and discuss the diversity-multiplexing tradeoff as 

well as the throughput-reliability tradeoff. From these results, we show that the proposed multi-stream cooperative communication 
scheme achieves lower outage probability compared to existing baseline schemes. 
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0> ! I. Introduction 

Cooperative communications for wireless systems has recently attracted enormous attention. By utilizing cooperation among 
different users, spatial diversity can be created and this is referred to as cooperative diversity |[T|-||3l. In HI, the authors 
considered the case where multiple relay nodes are available to assist the communication between the source and the destination 
^ ■ nodes using the decode-and-forward (DF) protocol, and they showed that a diversity gain of M(l — 2r) can be achieved with 
M relay nodes and a multiplexing gain of r. However, the advantages of utilizing more relay nodes is coupled with the 
consumption of additional system resources and power, so it is impractical (or even infeasible) to activate many relay nodes in 
resource- or power-constrained systems. As a result, various relay selection protocols have been considered in the literature. For 
■ example, in |l4l-||7l the authors considered cooperative MIMO systems with a single-antenna destination node and the selection 
of one relay node in the cooperative phase. In [4j the authors proposed an opportunistic relaying protocol and showed that this 
\^ ' scheme can achieve the same diversity-multiplexing tradeoff (DMT) as systems that activate all the relay nodes to perform 
distributed space-time coding. In IS) the authors applied fountain code to facilitate exploiting spatial diversity. In JS), Q the 
authors proposed a dynamic decode-and-forward (DDF) protocol which allows the selected relay node to start transmitting as 
soon as it successfully decodes the source message. On the other hand, there are a number of works lISl- lfTOl that studied 
the capacity bounds and asymptotic performance (e.g. DMT relation) for single-stream cooperative MIMO systems with a 
single-antenna destination node. In [Sj the authors derived the DMT relation with multiple full-duplex relays and showed that 
the DMT relation is the same as the DMT upper bound for point-to-point MISO channels. In ||9l the authors studied the 
network scaling law based on the amplify-and-forward (AF) relaying protocol. In all the above works, the destination node 
is assumed to have single receive antenna and hence, only one data stream is involved in the cooperative phase. When the 
destination node has multiple receive antennas, the system could support multiple data streams in the cooperative phase and 
this could lead to a higher spectral efficiency. 

In this paper, we design a relay selection scheme for multi-stream cooperative systems and analyze the resultant system 
performance. In order to effectively implement multi-stream cooperative systems, there are several technical challenges that 
require further investigations. 

. How to select multiple relay nodes to support multi-stream cooperation in the cooperative phase? Most of the existing 
relay selection schemes are designed with respect to having a single-antenna destination node and supporting cooperative 
spatial diversity. For example, in jlH the authors considered a rateless-coded system and proposed to select relay nodes for 
supporting cooperative spatial diversity based on the criterion of maximizing the received SNR. In order to support multi- 
stream cooperation, the relay selection metric should represent the holistic channel condition between all the selected relay 
nodes and the destination node, but this property cannot be addressed by the existing relay selection schemes. 
< How much additional benefit can multi-stream cooperation achieve? The spectral efficiency of the cooperative phase 
can be substantially increased with multi-relay multi-stream cooperation compared to conventional schemes. However, the 
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system performance may be bottlenecked by the source-relay links. Therefore, we characterize the advantage of multi-stream 
cooperation in terms of the end-to-end performance gain over conventional schemes that are based on cooperative spatial 
diversity. 

We propose a multi-stream cooperative relay protocol (DF-MSC-opt) for cooperative systems with a multi-antenna destination 
node. We consider optimized node selection in which a set of relay nodes is selected for multi-stream cooperation. Based 
on the optimal relay selection criterion, we compare the outage capacity, the diversity-multiplexing tradeoff (DMT), and the 
throughput-reliability tradeoff (TRT) of the proposed DF-MSC-opt scheme against traditional reference baselines. 

Notation: In the sequel, we adopt the following notations. C^^^^ denotes the set of complex M x TV matrices; Z denotes 
the set of integers; upper and lower case letters denote matrices and vectors, respectively; {■)'^ denotes matrix transpose; Tr(-) 
denotes matrix trace; diag(a::i, . . . , X]^) is a diagonal matrix with entries xi, . . . , x^, !(•) denotes the indicator function; Pr(X) 
denotes the probability of event X; \n denotes the N x N identity matrix; [•] and [-J denote the ceiling and floor operations, 
respectively; = max(-,0); =, < and > denote exponential equality and inequalities where, for example, f{p) ^ if 
limp^oo ^°\og(p)'' = ^' ^y['] denotes expectation over Y; and F{x,k) denotes the cumulative distribution function (CDF) 
for value x and degrees-of-freedom k. 



II. Relay Channel Model 

We consider a system consisting of a single-antenna source node, M half-duplex single-antenna relay nodes, and a destination 
node with antennas. For notational convenience, we denote the source node as the 0* node and the M relay nodes as 
the {1,2,..., M}-th node. We focus on block fading channels such that the channel coefficients for all links remain constant 
throughout the transmission of a source message. 

We divide the transmission of a source message that requires N channel uses into two phases, namely the listening phase 
and cooperative phase. In the listening phase, all the relay nodes listen to the signals transmitted by the source node until K out 
of the M relay nodes can decode the source messag^H In the cooperative phase, the destination node chooses Nr nodes from 
amongst the source node and the successfully decoding relay nodes to transmit multiple data streams to the destination node. 
Specifically, let x = [.t(1), x(2), . . . , x{N)Y' 6 C^^^ denote the signals transmitted by the source node over N channel uses. 
Similarly, let x.^ = [x„i{l), Xm(2), . . . , XmiN)]^ G C^^^, to = 1, ... , M, denote the signals transmitted by the m* relay 
node over channel uses. The signals received by the to*'' relay node is given by y„i = [j/,„(l), ym(2), . . . , ?/m(A^)]-^ £ C^^^, 
TO = 1, . . . , M, where 

Vrain) = hsR,mX{n) + Zyn{n), (1) 

hsR.rn is the fading channel coefficient between the source node and m*'' relay node, HsR = [hsR.i, hsR,2, ■ ■ ■ , hsR.Ai]'^ G 
(j^A/xi ^ vector containing the channel coefficients between the source and the M relay nodes, and Zm{n) is the additive noise 
with power normalized to unity. Each relay node attempts to decode the source message with each received signal observation 
until it can successfully decode the message. The listening phase ends and the cooperative phase begins after K relay nodes 
successfully decodes the source message. In the cooperative phase, let Xk{n), k — 1,...,K, denote the signal relayed by 
the fc*'' successfully decoding relay node and the aggregate signal transmitted in the cooperative phase can be expressed as 
^oin) = [x{n), xi{n), X2in), . . . £ £^{K+i)xi^ Accordingly, the received signals at the destination node are given 

by Y = [y(l),y(2),...,y(iV)]^ G C^^^-, where 



y(") 



HsD^in) + z{n) for the listening phase 

Hi3(X')VxD(n) + z(n) for the cooperative phase 

tlsD G C^'^^ represents the fading channel coefficients between the source and the destination nodes, T> is the set 
representing the K successfully decoding relay nodes, li]j{'D) = [HsoiilRD,!, '^rd.2, ■ ■ ■ , ^rd.k] G C^'^^'^+^-' represents 
the aggregate channel in the cooperative phase with Hj^xi.fc G C'^''^^ being the fading channel coefficients between fc*'* 
successfully decoding relay and the destination node, z{n) e c^rxi additive noise with power normalized to unity, and 

V is the node selection matrix. The selection matrix is defined as V = diag(wo, wi, W2, . . . , vk), where Vk = 1 if the fc* node 
(k G {0, . . . , K}) is selected to transmit in the cooperative phase and Vk = if the node is not selected. 
The following assumptions are made throughout the rest of the paper. 

Assumption 1 (Half-duplex relay model): The half-duplex relay nodes can either transmit or receive during a given time 
interval but not both. ■ 

Assumption 2 (Fading model): We assume block fading channels such that the channel coefficients UsR. and Hi3(P) 
remain unchanged within a fading block (i.e., N channel uses). Moreover, we assume the fading channel coefficients of 
the source-to-relay (S-R) links, relay-to-destination (R-D) links, and source to destination (S-D) links are independent and 
identically distributed (i.i.d.) complex symmetric random Gaussian variables with zero-mean and variance crl^, ctrd and (JgD, 
respectively. ■ 

'The relay system cannot enter the cooperative phase if less than K relay nodes can decode the source message within A' channel uses. 
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Assumption 3 (CSI model): Each relay node has perfect channel state information (CSI) of the Unk between the source node 
and itself. The destination node has perfect CSI of the S-D link and all R-D links. For notational convenience, we denote the 
aggregate CSI as H {Usr, 11d{V)). ■ 

Assumption 4 (Transmit power constraints): The transmit power of the source node is limited to ps- The transmit power of 
the k*^ relay node is limited to pk- ■ 

III. Problem Formulation 

In this section, we first present the encoding-decoding scheme and transmission protocol of the proposed DF-MSC-opt 
scheme. Based on that, we formulate the multi-relay selection problem as a combinatorial optimization problem. 



A. Encoding and Decoding Scheme 

The proposed multi-stream cooperation system is facilitated by random coding and maximum-likelihood (ML) decoding 
lfT2l . At the source node, an i?-bit message W, drawn from the index set {1, 2, . . . , 2^}, is encoded through an encoding 



function : {1,2,... 2^} — )■ X^. The encoding function at the source node can be characterized by a vector codebook 



C ^ {X^(1),X'"(2),...,X^^(2«)} G C 



rN 



The m* codeword of codebook C is defined as 



„(") 



(1) 



1^) 



4:^(1) 



(m) 



(Ni + l) 



(iVi + 1) 



„(") 



(m) 



(N) 



(N) 



, m 



(2) 



which consists of N vector symbols of dimension N,. x 1, and x)," (n) is the symbol to be transmitted by the k antenna 
during the n* channel use. The vector codebook C is known to all the M relay nodes (for decoding and re-encoding) as well as 
the destination node for decoding. At the receiver side (relay node or destination node), the receiver decodes the R-bit message 
based on the observations, the CSI, and a decoding function : (Y^ x H) — > {1, 2, . . . , 2^}. We assume ML detection in 
the decoding process. The detailed operation of the source node and the relay nodes in the listening and cooperative phases 
are elaborated in the next subsection. 



B. Transmission Protocol for Multi-Stream Cooperation 

The proposed transmission protocol is illustrated in Fig. [T] and the flow charts for the processing by the source, relay, and 
destination nodes are shown in Fig. |2] Specifically, the A^-symbol source message codeword (cf. (|2]i) is transmitted to the 
destination node over two phases; the listening phase that spans the first A^i channel uses and the cooperative phase that spans 
the remainder N2 = N — Ni channel uses. 

1 ) Listening phase: In the listening phase, the single-antenna source node transmits the first row0 of the message codeword 
(i.e. [x'(^\l) . . . x'(^\Ni)\ as per (|2]i) to the relay and destination nodes. Each relay node attempts to decode the source 
message with each received signal observation. Although the source node transmits only the first row of the source message 
codeword, the relay nodes can still detect the source message using standard random codebook and ML decoding argument. 
Effectively, we can visualize a virtual system with a multi-antenna source node as shown in Fig. |3] and the missing rows in 
the message codeword transmitted by the source node is equivalent to channel erasure in a virtual MISO source-relay channel. 
Once a relay node successfully decodes the source message, it sends an acknowledgement ACKbd to the destination node 
through a dedicated zero-delay error free feedback hnk. 

2) Control phase and signaling scheme: Without loss of generality, we assume that K relay nodes can successfully decode 
the source message with A^i received signal observations. Upon receiving the acknowledgement from the K successfully 
decoding relay nodes, the destination node enters the control phase and selects Nr nodes to participate in the multi-stream 
cooperation phase (cf. Section ITlI-Cl l. Specifically, the destination node indicates to the source and relay nodes the transition to 
the cooperative phase as well as the node selection decisions via an (M + l)-bit feedback pattern. The first bit of the feedback 
pattern is used to index the 0* node (the source node) and the last M bits are used to index the M relay nodes. The feedback 
pattern contains Nr bits that are set to 1 ; the m* bit of the feedback pattern is set to 1 if the corresponding node is selected 
to participate in the cooperative phase, whereas the bit is set to if the node is not selected. Note that the total number of 
feedback bits required by the proposed multi-stream cooperation scheme is K ACKnn plus one feedback pattern with 7\/ + 1 
bits, which is less than 2 bits per relay node. 

^The source node has a single transmit antenna and hence, could only transmit one row of the vector codeword {m) during the listening phase. 
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3) Cooperative phase: In the cooperative phase, the Nr selected nodes cooperate to transmit the (A^i + 1)* to the A^* 
columns of the source message codeword (cf. ^) to the destination node to assist it with decoding the source message. 
Specifically, for k — 1, . . . , Nr, the node corresponding to the k*'^ bit that is set to 1 in the feedback pattern transmits the fc* 



row of the message codeword (i.e. [x^^\Ni + 1) . ..x[!^>{N)] as per (|2])) to the destination node. 

To better illustrate the proposed transmission protocol, we show in Fig. [T] an example of the proposed system with a 
destination node with Nr = 2 antennas. Suppose the feedback pattern is given as follows: 

^n°oT 1^ A/ relay nodes ^| 

Feedback Pattern 






1 













1 



A/ + 1 bits 



This corresponds to selecting the first relay node (i?i) and the M relay node {Rm) to participate in the cooperative phase. In 



the listening phase, the source node transmits the first row of the codeword X (m) given by [x^;"\l)...xr'{Ni)]. In the 



cooperative phase, i?i and Rm transmit [x''"^\Ni + 1) . . .x'i^^N)] and [x^\Ni 
the effective transmitted codeword can be expressed as 



(m) 



(A^)], respectively. Therefore, 



Transmitted by source 
node in listening phase 



Transmitted by i?i 
in cooperative phase 




Transmitted by Rm 
in cooperative phase 



C. Problem Formulation for Multi-stream Cooperation 

We focus on node selection by the destination node for optimizing outage performance (cf. Section IIII-B2l i. where the 
destination node only has CSI of the S-D and R-D links 'Rd{V) (cf. Assumption O. The proposed Multi-stream Cooperation 
Scheme with Optimized node selection is called DF-MSC-opt and we define the outage event as follows; 

Definition 1 (Outage): Outage refers to the event when the total instantaneous mutual information between the source and 
the destination nodes is less than the target transmission rate R. Mathematically, outage can be expressed as 

-MSC-opt (H,V)<i?), (3) 

where H = {H.sr^'^d{T^)) is the aggregate channel realization and V is the node selection action given H£)(P). ■ 
The total instantaneous mutual information of the multi-stream cooperative system in (|3]l is given by the following theorem. 



Theorem 1: (Instantaneous Mutual Information of Multi-Stream Cooperative System) Given the aggregate channel realization 
H = (Hs/j, H£)(P)), the instantaneous mutual information /DF-MSC-opt(H, V) (bits/second/channel use) between the source 
and destination nodes for multi-stream cooperative system can be expressed as 

^DF^MSC-opt 

where A'^i and N2 are the numbers of channel uses of the listening phase and the cooperative phase, respectively, and F = 
diag{ps, pi, P2, ■ ■ ■ ,Pm) is the transmit power matrix of the source node and the M relay nodes. Note that iVi and N2 are 
random variables that depend on the realization of the S-R links Hsr- ■ 

The proof of Theorem [T] can be extended from |[T2l . ifTsl by applying random Gaussian codebook argument. Note that the first 
term in the mutual information in ^ corresponds to the contribution from the source transmissioiJ§ By virtue of the proposed 
DF-MSC-opt scheme, multiple data streams are cooperatively transmitted in the cooperative phase (unlike traditional schemes 
wherein only a single stream is transmitted). Hence, we can fully exploit the spatial channels created from the multi-antenna 
destination node and therefore achieving higher mutual information in the second term of (|4]i. 

By Definition [T] the average outage probability is given by 

VoutCV) - EH[X(/DF-MSC^opt(H, V) < R)] , 

which is a function of the node selection policy V. It follows that, given the channel realization H, the optimal node selection 
policy is given by V* = arg min l(/DF-MSC-opt(H, V) < i?), where the set of all feasible node selection actions is defined as 

n = {A e {0, i}(^^+i)x(Af+i)|A is diagonal and Tr(A) = Nr}. (5) 
Equivalently, for any transmission rate R, the optimal node selection policy is given by 

V* = argmax/DF-MSC-opt(H, V). (6) 



^The mutual information contributed by the source node can also be obtained from the mutual information of the virtual multi-antenna source model in 
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IV. Outage Analysis for the DF-MSC-opt Scheme 

In this section, we derive the asymptotic outage probability of the DF-MSC-opt scheme. For simplicity, we assume (jgjj = 
= crj-, and the power constraints of all the nodes are scaled up in the same ordeiQ, i.e., limps_i.oo 77 ^ 1 for 
m e {1, . . . , M}. The outage probability can be expressed as 

-MSC-opt (H,V*) <i?)] = EH,„[Pr(/DF -MSC-opt (H,^^) < i?|HsH)], (7) 

where in (a) the randomness of /DF-MSC-opt(H, V*) is introduced by the aggregate channel gain in the cooperative phase 
H£)(D) and the optimal node selection policy V*. Let Hni denote the collection of all realizations of the S-R links HsR 
such that K out of M relay nodes can successfully decode the source message, and the listening phase ends at the (A'^i)* 
channel use. We define the outage probability given IIsr and the duration of the listening phase A^i as PsR{iisB., Ni) = 
Pr(/DF-MSC-opt(H, V*) < RlHsR G Hni), which includes the following cases. 

« Case 1: The listening phase ends within N channel uses. In this case, Ni < N, N2 ^ N — Ni, and the outage probability 
Psr(H.sRi Ni) can be expressed as 

pC-i(Hsfl,iVi)=Pr[(^log(l + ps|HsDp) + ^5(H,V*)) <i?], (8) 

whereg(H,V*) = maxlogdet {lN,+psliDiV)VV"liD{V)") fovK > iV,., and ^(H, V*) = logdet {lN,+psliDiV)B.D{V)") 
for K < Nr. ^ 

• Case 2: The listening phase lasts for more than N channel uses. In this case, the outage event is only contributed by the 
direct transmission between the source and the destination nodes. The outage probability Psr^H-bRt^i) can be expressed 
as 

\B.sR, N{) = Pr( log(l + ps\iisD?) < R) . (9) 
Substituting ([8]l and (|9]l into O, the outage probability is given by 

Kut = Ell PsT ' (H5H, l)Pr{HsR eHi)+ J:i>n PsT 'i^SR, W^sr E Hi) (10) 
= E^Ii P§R ' (HsiJ, OPr(Hsfl e -Hi)+Pr{ log(l + ps|Hsi, < i?)Pr (Hsr e ^JiyNUi) , 

where Pr(H5fl & TLi) is the probability that K out of M relay nodes can successfully decode the message at the channel 
use. To evaluate the outage probability, we calculate each term in equation ( fTOl i as follows. First, it can be shown that 

Pr(Hsii e T^jvJ =Pr(Hsi? G UiyN^-iHi) - Pr (Hsr G Ui>n,Hi) = ^n,-i-<^n, (H) 

where 

- -exp(-^))--(cxp(-^))\ ^ . l,...,iV, (12) 

denotes the probability that less than K relay nodes can successfully decode the source message at the Z*'' channel use. Second, 
Pg^" '(Hsj7, 1) can be upper-boundecd as shown in the following lemma. 

Lemma 1: The outage probability for the DF-MSC-opt scheme given that K relay nodes can successfully decode the source 
message at the channel use can be upper bounded as 

Pl'j^^'iUsRj) <Pi{f{ai,H,V*) < R), l^l,...,N, (13) 

where ai = l/N, /(a,, H, V*) = ai log(l + ps\iisD\^) + (1 - «;) EfJi log2 (l + Ps<yh<^)) with Lt = min(iV,, K + 1) 
denoting the number of transmitted streams in the cooperative phase, and K{i), i = 1, . . . , Lt, are {Lt out of K + 1) ordered 

-distributed variables with 2Nr degrees of freedom. 
Proof: Refer to Appendix |A] for the proof. ■ 

Third, when fewer than K relay nodes can successfully decode the source message within N channel uses, the source 
node transmits to the destination node with the direct Unk only; it follows that Pr(H5fl; G U/>Ar'H;) = $Ar and Pr(log(l + 

PsIHsdP) < R) = F ( 2— ^;iVr ). Therefore, the outage probabiHty (cf. ([Toll) is given by 

-PLt = <PF[^;N,.)+j:l=i/Ni'^i-i-^i)P§T'iiisR,l) 

< <i>F (g^; A^.) + Ea,=i/^(1>;-i - $/)Pr (/(a/, H, V*) < R) 

= cl,F(g^;7V,) -ci>;,Pr(/(a;,H,V-)<i?), (14) 

*As such, there is no loss of generality to assume F = ps^M+i when studying high SNR analysis. 

'it is non-trivial to evaluate ' (lisR, exactly due to the dynamics of the optimal nodes selection. 
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where <i>J and ^J, denote the first order derivative of with respect to ai and evaluated at I = aiN and I' ~ ct'iN, respectively. 
Note that in the last step of ( fT4b we apply the Mean Value Theorem lfT4l which guarantees the existence of the point a[ £ (0, 1). 
Moreover, the outage probability upper bound of the DF-MSC-opt scheme is given by the following theorem. 
Theorem 2: For sufficiently large A^, the outage probability upper bound of the DF-MSC-opt scheme is given by 

A'+l 

rtut<^F[^-Nr]-^[,F{^-^-Nr]F{^-^^^-Nr] . (15) 



Direct Transmission Relay-Assisted Transmission 

Proof: Refer to Appendix |B] for the proof. ■ 
For systems without relays, the outage probabiUty is given by Vout.,pp — F (^j^-j^^ ; Nr^ . Note that the DF-MSC-opt 
scheme can achieve lower outage probability by taking advantage of multi-stream cooperative transmission. We can interpret 
— as the cooperative level, which quantifies the probability that the relay nodes can assist the direct transmission. For 
example, — $J/ =1 implies that the DF-MSC-opt scheme can be fully utilized, whereas =0 implies that the source node 
transmits to the destination node with the direct link only. As per (fT4l i. the outage probability Vout ^ decreasing function 
with respect tc|l— whereas — is a decreasing function with respect to the strength of the S-R links cr|j^. As shown in 
Fig. in when the strength of the S-R links increases from 30dB to 40dB, the cooperative level increases and hence the outage 
probability decreases. On the other hand, as shown in (flST l and as illustrated in Fig. |5] the outage probability decreases with 
increasing number of receive antennas at the destination node. 



V. DMT AND TRT Analyses for the DF-MSC-opt Scheme 

In this section, we focus on characterizing the performance of the DF-MSC-opt scheme in the high SNR regime. Specifically, 
we perform DMT and TRT analyses based on the outage probability V^^f as shown in ( fT4l i. There are two fundamental 
reasons for the performance advantage of the proposed DF-MSC-opt scheme, namely multi-stream cooperation and optimal 
node selection V* (cf. (|6]l). To illustrate the contribution of the first factor, we shall compare the performance of DF-MSC-opt 
with the following baselines: 

< DF-SDiv (Baseline 1): The DF relay protocol for cooperative spatial diversity. The listening phase and the cooperative phase 
have fixed durations, i.e., each phase consists of N/2 channel uses. The source node and all the successfully decoding relay 
nodes cooperate using distributed space-time coding. At the destination node, maximum ratio combining (MRC) is used to 
combine the observations from the different receive antennas. 

< AF-SDiv (Baseline 2): The AF relay protocol for cooperative spatial diversity. All the relay nodes transmit a scaled version 
of their soft observations in the cooperative phase. At the destination node, MRC is used to combine the observations from 
different receive antennas. 

< DDF (Baseline 3): The dynamic DF protocol /|6^. Once a relay node successfully decodes the source message, it immediately 
joins the transmission using distributed space-time coding. At the destination node, MRC is used to combine the observations 
from different receive antennas. 

On the other hand, to illustrate the contribution of optimal node selection, we shall compare the performance of the DF-MSC- 
opt scheme with the DF-MSC-rand scheme (Baseline 4), which corresponds to a similar multi-stream cooperation scheme but 
with randomized node selection V randomly generated from the node selection space O (cf. (|5]l). Table U summarizes the major 
differences among the DF-MSC-opt scheme and the baseline schemes. We illustrate in Fig.|6]the outage capacity versus SNR 
(with Vout = 0.01) of the cooperative systems with Nr = 3, K = 3, M = 15. Note that the proposed DF-MSC-opt scheme 
can achieve a gain of more than 1 bit/channel use over the four baseline schemes. 



A. DMT Analysis 

In order to analyze the DMT relation of the DF-MSC-opt scheme, we first derive the relation between the outage probability 
and the multiplexing gain. In the outage probability expression (fT4] i. for a sufficiently large ps, the asymptotic expression of 
the first term is given by 



and the asymptotic expression of <i>J, is given by 

= Ef=o'(f)P5 ' <PS ' =PS ' (17) 



''For a sufficiently large SNR ps, F{ ; A^r ) F ( ^^^'^ y'' ~^ ■,Nr] is much smaller than F ( J'^~'^'2 ; A^r 



K+l 



TO APPEAR IN IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, 2010 



7 



where (j>i' = exp I — ^ ^J^.;, j . Moreover, we can express the term Pr (/(aj, H, V) < R) in ( fT4l i as 

Pr(/(a;,H, V) < i?) - Pr(pf OfJi P^-^^*^-^'^' < P^s) 

= Pr(a;(l - S)+ + EfJi(l - - < ^ //e^l'^, /3)d^d/3, (18) 

where —5 is the exponential order of cr|j), —pi is the exponential order of 7(1), B = {5, /3 : Q!;(1 — (5)+ + (l — aj) < 
r}, and p((5, /3) is the joint probability density function of 6 and /3. Substituting (fT6]l-(fr8]l into (fT4l i. the outage probabiUty can 
be expressed as 

VL. < p,(---+^+-^)(-^)%pr'"""^^^"^^'//,p(J,/3)dMA (19) 

and the DMT relation for the DF-MSC-opt scheme is given by the following theorem. 
Theorem 3: The DMT relation for the DF-MSC-opt scheme can be expressed as 

d{r, i^)DF-MSC-opt = min((ii, ^2 + da), (20) 

where di = (M -K + 1 + iV^)(l - r)+, da = (M -K + 1)(^)+, and 

+ (iV^ - 0<r<l/2 
d-^ ■ ( N,.+{Nr-e){K+i-9) {N^-e)(K+i-e)\ ^ < e+i 

«3 - S nmn^_(jv^^^_2g)(_^_e), J S ? < (/I) 

with 6* = 1, 2, . . . , Lt and I/T = mm{Nr, K + 1). 

Proof: Refer to Appendix |C] for the proof. ■ 
We could further optimize the parameter K in the DF-MSC-opt scheme and the resulting DMT relation is given by 

'^(^)DF-MSC-opt = , max d{r, A')DF-MSC-opt =max {d{r, lK*\ )DF-MSC-opt, d{r, \K*] )DF-MSC-opt) , 



where 



(M+l){2-4r+r^) + Nr.{2-3r+r^) 



< r < 1/2 



I (A/+l)(l-r)^-Af^(3r-r^-(l-r)2e)-e(36l(l-r)-l-r) fl , ^ 6 + 1 

{ Nr(l-r)-r+{l-ry^ 6 + 1 ^ ' — 8+2 

is the optimal number of successfully decoding relay nodes to wait for in the listening phase. 

In the following, we show that the DF-MSC-opt scheme achieves superior DMT performance than the traditional cooperative 
diversity schemes. Specifically, as per |T5|, the AF-SDiv and DF-SDiv schemes have identical DMT relations given by 
c?(r)AF/DF^SDiv = - 2r)+ + Nr{l - r)+ for < r < 1. For the DDF protocol, the DMT relation is given by the 

following lemma. 

Lemma 2: The DMT relation for the DDF protocol with Nr receive antennas at the destination node can be expressed as 

r {M + Nr){l-r) 0<r<i 
d(r)DDF-<^ Nr + M{^)+ i<r<i (23) 
I NrC-^) \<r<l 

Proof: Refer to Appendix |D] for the proof. ■ 
In Fig. |7] we compare the DMT relations for the DF-MSC-opt scheme and the baseline schemes. For a system with M = 15 
relay nodes and a destination node with iV^ = 3 antennas. Note that since the DF-MSC-opt scheme (as well as the DF-MSC- 
rand scheme) exploits multi-stream cooperation, it can achieve high diversity gain than the traditional cooperative diversity 
schemes (i.e. AF-SDiv, DF-SDiv, and DDF) that exploit single-stream cooperation. Moreover, since the DF-MSC-opt scheme 
optimizes the node selection policy and the number of decoding relay nodes K, it can achieve better diversity gain than the 
DF-MSC-rand scheme. 



B. TRT Analysis 

The DMT analysis alone cannot completely characterize the fundamental tradeoff relation in the high SNR regime. Specifi- 
cally, the multiplexing gain gives the asymptotic growth rate of the transmission rate R at high SNR p and is only applicable 
to scenarios where R scales linearly with log p. Hence, there are many unique transmission rates that correspond to the same 
multiplexing gain, and the DMT analysis only gives a first order comparison of the performance tradeoff at high SNR when 
we have different multiplexing gains. 
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In order to have a clearer picture on the tradeoff relations, in the following we study the SNR shift in the outage probabiUty 
Tout as we increase the transmission rate R by AR. We quantify the more detailed relations among the three parameters 
{R,\ogp,Vout{R,p)) by analyzing the TRlQ relation il6 |. 

Consider the outage probability expression for the DF-MSC-opt scheme ( fT4b . In order to facilitate studying the TRT, we 
choose the parameter K such that the outage events are dominated by relay-assisted transmissions, i.e. K > \K*^ (where K* 
is given by (l22b). The following theorem characterizes the asymptotic relationships among R, ps, and for the t > \ 
cas^ 

Theorem 4: The TRT relation for the DF-MSC-opt scheme under K > \K*~\ and r > i is given by 

'°g^°"'7;"°P-(-)^ = -gDF-MSC-opt(2) (24) 

where 

n{z) ^ {{R, ps)\z + 1 > >z} for zeZ,0<z<LT, (25) 

denotes the 2* operating region, coF-MSC-optC^) = A' + 1 + N,, - (2z + 1), and .gDF^MSC-opt(^) ^ (K + l)Nr - z{z + 1). Note 
that S(DF-MSC-opt(2) is defined as the rehabihty gain coefficient and tDF-MSC-opt(^) = .gDF-MSC-opt(^)/cDF-MSC-opt(z) is defined as 
the throughput gain coefficient. 

Proof: Refer to Appendix |E| for the proof. ■ 
By applying TheoremlH the SNR shift between two outage curves with a AR rate difference is 3AR/t{z) dB. Fig.|8]shows 
the outage curves corresponding to AR = 2 bits/channel use for Nr = 3, K = 3 and M = 15. This scenario corresponds to the 
region 7^(1) and the TRT relation for the DF-MSC-opt scheme can be expressed as (?DF-MSC-opt(l) = 10 and iDF-MSC-opt(l) = §■ 
As we can see from the simulation results, the SNR shift is 2.4 dB for 2 bits/channel use increase in the transmission rate, 
which matches with our analysis. 



VI. Conclusion 

In this paper, we proposed a multi-stream cooperative scheme (DF-MSC-opt) for multi-relay network. Optimal multi-relay 
selection is considered and we derived the associated outage capacity as well as the DMT and TRT relations. The proposed 
design has significant gains in both the outage capacity as well as the DMT relation due to (1) multi-stream transmissions in 
the cooperative phase and (2) optimized relay selection. 



Appendix A 
Proof of Lemma [T] 

To find the outage probability given I, we characterize the function of g(H, V*) under different conditions. 

• Condition 1 K < Nri Under this condition, we allow all the successfully decoding relays to transmit in the cooperative 
phase. Thus, the communication links can be regarded as a conditional MIMO link ifTTl and ^(H, V) = log dot (In^ + 
PsHd(2?)H£,(P)^) < J2fJi^ log2(l + PsCTofiii)), where K{i), i 1, . . . , A' + 1, are x^-distributed variables with 2Nr 
degrees of freedom. 

• Condition 2 K > Nr'. Under this condition, we select Nr nodes out of the source and the successfully decoding 
relay nodes from the node selection space $1 to participate in cooperative transmission. The analytical solution for the 
term g(H,V) is in general not trivial i fTSll . Since in the proposed DF-MSC-opt scheme we choose the best Nr out of 
A' + 1 transmit nodes, the upper bound can be obtained similar to IfTSl (6)]. Thus, the capacity bound with transmit node 
selection is .g(H, V) = maxvGfJ logdet psliD{V)VV"llD{V)") < log2(l + ps(jj)K{i)), where the K{i), 
i = 1,. . . ,Nr, are ordered x^-distributed variables with 2Nr degrees of freedom. 

Combining the two conditions above, we obtain the general form as shown in Lemma [T] 



Appendix B 
Proof of Theorem[2] 

From (O, the only unknown part is given by 

Pr(/(a;, H, V) < i?) = Pr[a; log(l + psIHsd?) + (1 - a\) Efji log2(l + Ps<yh<i)) < R] , 

'TRT analysis allows for investigating more general scenarios where the transmission rate R does not scale linearly with log p and helps to study the SNR 
gain of the outage probability vs SNR curve when we increase the transmission rate _R by AR. 

*For the f < 5 case, the MIMO channel formed by the multiple relay nodes to the destination node can only operate with a multiplexing gain of 1 and 
the TRT relation is trivial. 
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which can be relaxed as 



Pr[a\ log(l + psinsD^ + (1 - a[) EfJi log2(l + Ps^|,«(*)) < R] 
< Pr(«; log(l + psIHsDp) < i?)Pr((l - a[) log2(l + ps'Jl<^)) < R) 

R 

= F{^;Nr)Pr{ Efji log. (1 + PsalKi^)) < ^) . 
Moreover, the following relation holds for a[ > 

PrfEfJi log2 (1 + Ps^h^i^)) <T^)< Pr(log2 (1 + ^5^2,^(1)) < j^' 



where in the last step we applied the results of order statistics [[T9l . Substituting ( |26] | into ( fT4] i. we have Theorem |2l 

Appendix C 
Proof of Theorem[3] 

The asymptotic expression of the outage probability "Po^^ is given bv (fT9]l with B = {5,(3 : <5)+ + (l— aj) Ei(l~A)''' < 

r}. 

From the first term in (fT9] l. we can obtain di = {M — K + I + Nr){l — 

Since the source node is equipped with a single antenna, the maximum multiplexing gains for the S-R links as well as 
the S-D and R-D links is always less than 1. Specifically, the multiplexing gain in the listening phase is < 1 and the 
multiplexing gain in the cooperative phase is < 1 or equivalently, aj > max(r, 1 — r). Thus, we can evaluate the second 
term through = - A' + 1)(1 - j^) for r < 1/2 and otherwise. 

We evaluate the third term under the following cases. When r < 1/2, > 1 — r and i3 C {(5, /? : ?'((1 — <5)+ + Ei(l^ A)^) < 
r}. Equivalently, 

JJ^piS,(3)d6d(3 < //(^_,^+^^^(^_^^)+<^p(5,/3)d<5d/3 

— Pg ~ Ps ■ v^'/ 

As a result, da = {Nr - l)K + Nr for r < 1/2. When 1/2 < r < 1, > r and the corresponding B = {5,13 : 
r(l - S)+ + (1 - r) - l3r)^r < r}. Equivalently, 

JJ^piS,P)d6dP < //(i_,)++i^5^_(i_^^)+<iP(<5,/3)d<5d/3 

— Ps \-^°> 

for fz^ e {1, ■ . ■ , Lt}- Thus, we can write the general form as ^ info<,.'<i {Ny + (Nr - e){K + I - 0) ~ {Nr + K - 
29){j^ - 6)) for 6* = {1, 2, ... , Lt}- Since ds is hnear with respect to r', we can write d^ = min {Nr + {N,, - 6){K + 1 - 
61) - [Nr + K- 2e){-^ - e), Nr0i^) + {Nr - e){K + 1-6*)}. Moreover, when j^, > Lt, we have - Lt = and 

d^^NrLTi^). 

Combining the three terms above, we have Theorem [3] 

Appendix D 
Proof of Lemma|2] 

The main steps of the proof are based on the work by JS). Denote gk.j to be the channel coefficient between the j*'' node and 
the fc*'' receive antenna at the destination node, and the received signal at the fc*'' antenna is given by yk = gk,jXj+Zk- Following 
the same argument in the proof of Theorem 6 in ||6|, we can find that destinations ML error probability assuming error-free 
decoding at all of the relays, provides a lower bound on the diversity gain achieved by the protocol and the corresponding PEP 

for fc*'* antenna is upper bounded by < HjH^ ^ + ( Ei=i l5fe,iP)p > where rij is the number of symbol intervals in 
the codeword during which a total of j nodes are transmitting, so that, the total number of codeword length Ejfa^ = 
Hence, the PEP for the destination to received the information is upper bounded by P < Ofc^i R^- Note that, the upper bound 
is not tight since we do not consider the joint decoding among the receive antennas, but it is sufficient to derive the order-wise 

relation. By changing the order of k and j, we have P < 11^=1^ Ill^i + (Ei=i l3fc,iP)p) ■ X]j=i number 
of symbol intervals that relay p has to wait, before the mutual information between its received signal and the signals that the 
source and other relays transmit exceeds NR. Thus, we have Ej=i "-i — min|A^, iog{i+ \ h^+i iPp) }' where /ip+i,i denotes 
the channel condition between the relay node p and the source node. 
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Define Vk.j and Uj^i as the exponential orders of gk.i and hj^i. We have the following relation, 
where Vi = ■ Choose the rate R = r logp and the PEP bound for rate R is given by 

Hence, the set of channel realizations that satisfy {X]ji=i ^ (1 ^ min{ui, . . . , Vj})^ < r results in the outage event. 
Define iX, ~ min{wi, . . . , Uj} and we can separate the discussion in the following three cases. (Naturally, we have v\ > 

V2 > . .. > v~m). 

• Case 1; ?7i < 1. In this case, we can follow the discussion from (62) to (71) in ||6l and show that the diversity order d is 

r (iV, + M)(l-r), ^>r>0 

d>{ Nr + M{^), \>r> 



N,.+M ■ 

Case 2: u^v/ > 1. It's trivial and d = M. 



• Case 3: Vi> \> Ui+i. Follow the same discussion from (72) to (82) in ||6l , we conclude that d > 
Combining the above cases, we have Lemma |2] 



Appendix E 
Proof of Theorem|4] 

Given r > i and A' > \K*], ^ can be simplified as 7^*„t< - Pr(/(aJ, H, V*) < R). Substituting the expression of 
/(a;,H,V), we have 

V*,, < Pr(a; log(l + ps\iisD\') + (1 - a[) ^fji log2(l + Ps^!,7(«)) < R) 

< Pr((l-r)Efjilog2(l + P5fT^7«)<i?). (29) 
Based on the above results and following the same lines as the proof of |16, Theorem 2], we have the result of Theorem 2] 
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TABLE I 

Comparisons among the DF-MSC-opt scheme and the baseline schemes. 



Relay protocol 


AF-SDiv 


DF-SDiv 


DDF 


DF-MSC-rand 


DF-MSC-opt 


No. of transmitted streams in 
the cooperative phase 


1 


1 


1 


Nr 


Nr 


No. of relay nodes in the coop- 
erative phase 


M 


depends on 
S-R links 


depends on 
S-R links 


Nr 


Nr chosen 
from K 


Receiver structure 


MRC 


MRC 


MRC 


ML 


ML 
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f^eniu CDD^ native pLiase. Since tiie 
jed 1 ' is Hie secoiid bit-l' in tlie 
f&eclbnck patteiit. Km wiIL tiaLi^iiiit tlie 
^coiidi ow of tlie codewor<liu the 
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results ijqtificatioi 



one bit for SDurcE M fcrits, one for each rela^r 



Fig. 1. Illustration of timing diagram and an example of feedback pattern design. In this example, two relay nodes Ri and Rm are selected to transmit 
in the cooperative phase. The selected relay nodes will re-encode the same codeword (codeword X^(m) in this example) selected from the same common 
vector codebook C and the R-hit message received in the listening phase. Ri will send out the "first" row of codeword X^(m) and /2jv/ will send out the 
"second" row. 
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Fig. 2. Flow chart of the protocols at the source node, the relay node and the destination node. 




Fig. 3. The single-antenna source-relay channel is equivalent to a multi-antenna virtual MISO channel. 
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SNR (dB) 



Fig. 4. Outage probability vs. SNR of the DF-MSC-opt sclieme for different o-| „ under A^r = 3, X = 6 and M = 15. 




Fig. 5. Outage probability vs. SNR of the DF-MSC-opt scheme for different Nr under o-|^ = lOdB, K = G and M = 15. 




Fig. 6. Outage capacity comparison of different schemes for Nr = 3, K = 3 and A-I = 15. The channel variances of the S-R, R-D, S-D links are normalized 
to unity. 
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Diversity Multiplexing Tradeoff for Different Relay Protocols 




rvlultiplexing Gain r 



Fig. 7. Diversity-multiplexing tradeoff comparison of different relay protocols. 



Outage Probability vs. Total Transmit Power 




Total Transmit Power (dB) 



Fig. 8. Monte-Carlo simulation results for outage curves corresponding to AR = 2 bits/channel use for A^r = 3, K = 3 and M = 15 case with normalized 
S-R, R-D, S-D channel variances. 
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